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dr ABSTRACT 



o 



High resolution optical and ultraviolet spectra of two B-type post- Asymptotic Giant 
^ . Branch (post-AGB) stars in globular clusters, Barnard 29 in M13 and ROA 5701 in 

^ ' uj Cen, have been analysed using model atmosphere techniques. The optical spectra 

! have been obtained with FEROS on the ESO 2.2-m telescope and the 2d-Coude spec- 

trograph on the 2.7-m McDonald telescope, while the ultraviolet observations are from 
the GHRS on the HST. Abundances of light elements (C, N, O, Mg, Al and S) plus 
Fe have been determined from the optical spectra, while the ultraviolet data provide 
additional Fe abundance estimates from Fein absorption lines in the 1875-1900 A 
| wavelength region. A general metal underabundance relative to young B-type stars 

is found for both Barnard 29 and ROA 5701. These results are consistent with the 
£N) ' metallicities of the respective clusters, as well as with previous studies of the objects. 

The derived abundance patterns suggest that the stars have not undergone a gas-dust 
separation, contrary to previous suggestions, although they may have evolved from 
' the AGB before the onset of the third dredge-up. However, the Fe abundances de- 

, rived from the HST spectra are lower than those expected from the metallicities of 

the respective clusters, by 0.5 dex for Barnard 29 and 0.8 dex for ROA 5701. A sim- 
ilar systematic underabundance is also found for other B-type stars in environments 
, of known metallicity, such as the Magellanic Clouds. These results indicate that the 

■ Fe in ultraviolet lines may yield abundance values which are systematically too low 

by typically 0.6 dex and hence such estimates should be treated with caution. 
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1 INTRODUCTION 

Post- Asymptotic Giant Branch (post-AGB) stars are short- 
lived objects (typically 10 4 years), of initial mass between 
0.8 - 8 Mq. They have evolved from the AGB, under- 
gone severe mass loss , and will become planetary nebulae 
l|lben fe Renzinil Il983l ). T his evo l utiona ry st age has been 
review ed by, for example, iKwokl (| 1993f ) and IVan Winckell 
(2003), and is important as it provides an insight into the 
nucleosynthesis and mixing processes during the later stages 
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of stellar evolution. Most post-AGB stars have been identi- 
fied via their infrared excess, arising from a circumstellar 
dust shell. This shell is most readily detected at the earliest 
stages of post-AGB evolution, when it is relatively hot and 
close to the star. As a result, the majority of post-AGBs 
found to date have just left the AGB, and are of A, F or G 
spectral type (|Oudmaiierlll996l ). 

A number of abundance studies have been performed 
for post-AGB stars , with the majority being metal-poor 
i|Van Winckell I2003I ). However, there are many features 
of post-AGB evolution which are not yet fully under- 
stood. For example, several stars show abundance pat- 
terns reminiscent of interstellar gas, with Si, Ca etc fol- 
lowing the very low Fe abundance, while CNO and S 
are near solar (|Van Winckel. Waelkens fe Waters] Il995l ). 
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IWaters. Trams fe Waelkensl (|l992h proposed a mechanism 
for the formation of such a 'cleaned-up' photosphere in a 
post-AGB star, where grain formation occurs in the circum- 
stellar shell produced in the AGB phase of the star, and 
subsequent slow accretion of the gas from this material re- 
sults in the peculiar abundance pattern. This is most likely 
to happen during mass transfer between the components of 
a binary sys tem, or accretion tri g gered by a binary compan- 
ion. Indeed, IVan Winckel et "all (| 1995h showed that all the 
extremely Fe deficient post-AGB stars known at that time 
are binaries, providing support for this hypothesis. 

Although most of the known post-AGBs are cool ob- 
jects, several B-type stars originally classified as being Pop- 
ulation i are now believed to be hot post-AGBs, on the ba- 
sis of their derived atmospher ic parameters and chemical 
compositions (see, for example, iRvans et al.ll2003l and refer- 
ences therein). However, although the B-type stars should 
be evolving from the cooler post-AGBs, there are significant 
differences. The most striking is the large C underabundance 
found for several hot objects, which m ay imply that thes e 
left the AGB before the third dredge-up (fkvans et al.ll2003h . 

These apparent abundance differences between the hot 
and cool post-AGB candidates indicate that further investi- 
gations are required, so that this evolutionary stage may be 
more fully understood. As a first step, it is clearly important 
to study hot post-AGB objects of known initial metallic- 
ity, rather than field stars of unknown origin. Fortunately, a 
sample of these is available: the well-known UV-Bright stars 
found in some globular clusters, whose post-AGB nature has 
been unambiguously established fr om their position in t he 
cluster colour-magnitude diagram (|Landsman et all [2000'). 

A review of UV-B r ight s tars in globular clusters can 
be found in iMoehlerl l|200ll) . whi le searches for globu- 
lar cluster post-AGB stars include IZinn. Newell fe Gibson! 
1 19721 ) . iHarris. Nemec fe Hesserl l| 19831 ) and lde Boeij (|l985l v 
Of the B-type post-AGB stars studied to date, only 
three have well determ ined iron abundanc e estimates, 
namely ZNG-1 in M 10 dMoonev et al.l |2004), Barna rd 29 



in M13 (|Moehler et~al"1ll998l; iDixon fe Hurwitd 119981) and 



ROA5 701 in wCen (|Moehler et al.l Il998l ). iMoonev et all 



|2004 ) examined optical spectra of ZNG-1, finding an 
iron abundance larger than the metallicity of the clus- 
ter, M10. Barnard 29 and ROA5701 have been studied 
several times, primarily due to their relative brightness, 
using both optical and ult r aviolet spectr a (see, for ex- 



1984:; de Boer 


1983; Conlon. Dufton & Keenan 1994 


Adelman et al.l 


1994 Dixon & Hurwitzlll998l; iMoehler et all 


1998; Thompson et al. 2006). However the iron abundances 



have previous ly been estima t ed us ing only the UV spectra. 
In particular, IMoehler et al.l (|l998T ) have analysed HST ob- 
servations of these stars and derived Fe abundances which 
are lower than the cluster metallicities, which they ex- 
plained by the grain formation process discussed above. If 
Barnard 29 and ROA5701 have such 'cleaned-up' photo- 
spheres, then one would expect them to show similar deple- 
tions in Mg and Si, which are also subject to grain forma- 
tion. However, previous studies of Barn ard 29 and ROA 5701 
indicate no depletions in Mg nor S i llCqnlon et al l Il994 
IMoehler et al.lll998l ; Prhompson et aill2006h . 

Clearly, the abundances of Barnard 29 and ROA 5701 
need revisiting, in particular those for Fe. In this paper, high 



resolution optical and ultraviolet spectra of Barnard 29 and 
ROA 5701 are analysed, using model atmosphere techniques 
to estimate atmospheric parameters and chemical composi- 
tions. The observations and data reduction are discussed in 
Section 2, while Section 3 details the analysis methods and 
results. In Section 4 we present a discussion of our results, 
and in particular we assess the accuracy of Fe abundance 
determinations for B-type stars using ultraviolet spectra. 



2 OBSERVATIONS AND DATA REDUCTION 

Optical spectra of Barnard 29 and ROA 5701 were obtained 
during two observing runs, the former with the 2.7-m Har- 
lan J. Smith Telescope (McDonald Observatory) on three 
nights between 2005 April 30 and May 5, and the latter 
with the MPG/ESO 2.2-m Telescope (La Silla, Chile) on 
2006 April 10. Ultraviolet spectra of both stars were taken 
from the Hubble Space Telescope (HST) archive. Details of 
the observations are provided below. 



2.1 Optical spectra 

2.1.1 Barnard 29 

Barnard 29 was observed with the 2d-Coude spectrograph 
(|Tull et al.ll 19951 ) . using a 2.4 arcsec slit to achieve a spectral 
resolution (_R=A/AA) of ~35000 and a wavelength coverage 
of ~3700-10000 A. Due to cloud, data were taken only on 
three nights, resulting in 30 half hour exposures, or 15 hours 
on-source. A Th-Ar hollow cathode lamp was observed at 
the beginning and end of each night, with the centres of 
the lines stable to better than ~0.1 pixel. Test exposures 
were taken at the start of the run, to ensure th at optical 
ghosts ('the picket fence'; ICheng fe Ch ang 2000:) and bad 
columns did not overlap with wavelength regions of interest, 
in particular the Fein lines around 4420 A. Also, in order 
to ensure that the weak absorption line features detected 
were real, the central wavelength was changed, by ~0.4 A, 
after the first night and remained in the new position for 
the other two clear nights. 

The data were reduced using standard method s within 
the I mage Reduction and Analysis Facility (iraf; iBarned 
1993). This reduction included bias subtraction and flat- 
fielding, with stellar extr action and co smic-ray removal per- 
formed using doecslit (|Valdes|[l993l ) and scombine. The 
signal to noise ratio (S/N) per pixel in the co-added spec- 
trum was ~130 per pixel at ~4420 A (corresponding to ~180 
per resolution element), and this w as read into the Starlink 
package DIPSO l|Howarth et al„ 2004) for further analysis (see 
Section [2X3J. 



2.1.2 ROA 5701 

The optical spectra of ROA 5701 were obtained using the 
high-resolution Fibre-fed Extended Range Optical Spectro- 
graph (FEROS; iKaufer et al.lll999l ). with the EV 2Kx4K 
CCD detector and a 79 lines mm -1 echelle grating. In total, 
ROA 5701 was observed for 9 hours, resulting in a spectral 
resolution of R ~46000 and a wavelength coverage of ~3560- 
9200 A. 

The data were reduced on-line, using the standard 
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FEROS Data Reduction System (DRS), implemented with 
ESO-MIDAS software. This reduction includes the extrac- 
tion and de-blazing of the echelle orders, wavelength calibra- 
tion and merging of the orders. The reduced spectral expo- 
sures were combined into a single spectrum using SCOMBINE 
within iraf, and this was imported into DIPSO for further 
analysis (see Section |2 . 1 . 3 p - An average S/N ratio of ~85 
per pixel at ~4420 A was obtained, corresponding to ~120 
per resolution element. 



2.1.3 Equivalent width measurements 

Initially, the continua of the optical spectra for both 
Barnard 29 and ROA 5701 were normalised using low-order 
polynomials fitted to spectral regions free from absorption 
features. The line fitting program elf within DIPSO was then 
employed to measure absorption line equi valent widths, by 
fittin g Gaussian profiles to the lines (see iThompson et al] 
120061 for further details). Balmer lines and diffuse He I line 
equivalent widths were not measured, as instead theoretical 
profiles were fitted to the normalized data (see, for example, 
Fig. P. 



2.2 Ultraviolet spectra 

The ultraviolet spectra of Barnard 29 and ROA 5701 were 
obtained directly from the HST Multimission Archive 
( |http: / /archive.stsci.edu] z3EC0204T and z3EC0104M re- 
spectively). They were observed with the Goddard High 
Resolution Spectrograph (GHRS) using the G200M grat- 
ing (R = 25000) and the large science aperture (LSR), 
giving a wavelength coverage of ^18 60-1906 A. FP split s 
were not used for these observations l|Moehler et al.l l~1998). 
Barnard 29 was observed on 1996 November 30 for 4570 s, 
and ROA 5701 on 1996 August 3 for 4711 s. 

Following the standard pipeline reduction, the spectra 
were co-aligned within iraf, using the STSDAS package (Ver- 
sion 3.4). The combined spectrum for each star was then 
read into DIPSO and cross-correlated with theoretical spectra 
to obtain the stellar LSR radial velocities of -228 and +258 
km s _1 for Barnard 29 and ROA 5701, respectively. These 
are consistent w ith the radial velocities of M 13 (-227. 2 km 
s- MHarrisll 199fj ) ando;Cen (+240 km s-\ lNorris]| 19741). An 
inter stellar reddening correction, of E (B-V) = 0.02 (jHarrisI 
1 19961 ) for M13 and E(B-V) = 0.11 (|Cacciari et al.l Il984j ) 
for uj Cen, was applied to the appropriate spectrum. The 
spectra were then binned, with a wavelength interval of 0.05 
A, resulting in a wavelength coverage of 1868-1906 A for 
Barnard 29 and 1865-1903 A for ROA 5701. 



3 DATA ANALYSIS AND RESULTS 

3.1 Model atmosphere calculations 

Non-Local Thermodynamic Equilibrium (non-LTE) model 
atmosphere grids, calculated using tlusty and synspec 
dHubenvll988l ; lHubeny fc Lanzl 19951 ; iHubenv. Heap fc Land 
1998), have been used to derive atmospheric parameters 
and chemical abund ances. A summary of th e procedures in- 
volved are given in IThompson et al" (|2006h . while a more 
detailed discussion of the grids and methods may be found 
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Figure 1. Observed and theoretical line profiles for H7 at 4341 A 
in the spectra of Barnard 29 (upper panel) and ROA 5701 (lower 
panel). The theoretical profiles have been generated for surface 
gravities of 2.75, 3.0 and 3.25 dex (upper, middle and lower 
smooth curves, respectively) for Barnard 29, implying a surface 
gravity of 2.95 dex. For ROA 5701, results have been calculated 
for surface gravities of 3.0, 3.25 and 3.5 dex (upper, middle and 
lower smooth curves, respectively), implying log g = 3.3 dex. 



in lRvans et~aT1 (120031 ) and lDufton et all (|2005l ). and also at 
|http:/ /star, pst.q ub.ac.uk 

Briefly, four grids have been produced with metallicities 
corresponding to [Fe/H] = 7.5 dex for the Galaxy, 7.2 dex 
and 6.8 dex for the Large (LMC) and Small (SMC) Mag- 
ellanic Clouds, respectively, and 6.4 dex for low metallic- 
ity regimes. For each metallicity grid, approximately 3000 
models were calculated for a range of effective tempera- 
tures from 12000 to 35000 K, logarithmic gravities (log g 
= 4.5 dex down to the Eddington limit) and microturbu- 
lences of to 30 km s _1 . The iron abundance was fixed at 
each grid, while the light element abundances (C, N, O, Mg, 
Si, S) were allowed to vary from -0.8 to +0.8 dex around 
their base (metallicity) value. Iron is assumed to dominate 
the metal line blanketing. However, the exclusion of nickel 
in our models should not be a significant source of error 
in the atmospheric parameters and chemical composition 
l|Hubenv et al.| [l998') . The atmospheric structure is therefore 
defined by the iron abundance (metallicity), effective tern- 
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Table 1. Adopted atmospheric parameters for Barnard 29 and 
ROA 5701. 



Star 


Teff 


log g 


£ v sin? 




10 3 K 


dex 


km s _1 km s" 1 


Barnard 29 


20 ± 1 


2.95 ± 0.1 


4 ± 2 


ROA 5701 


25 ± 1 


3.3 ± 0.1 


1 ±\ 3 



perature, gravity and microturbulence. It is also assumed 
that light element abundances can vary with negligible affect 
on the atmospheric str ucture, as shown in tests performed 
bv lDufton etaH lj2005l ). 

Theoretical profiles and equivalent widths of light ele- 
ment transitions were calculated using the models. A GUI 
interface in IDL allows the user to access the theoretical 
equivalent widths and interpolate these to calculate equiva- 
lent widths and abundance estimates for ~200 metal lines, 
for any set of atmospheric parameters. Abundance incre- 
ments of 0.4 dex used in the grids were found to be suffi- 
ciently small that no significa nt errors were incu rred due to 
the interpolation procedures l)Rvans et al.l [2003). 



3.2 Atmospheric parameters 

The atmospheric parameters of each star were determined 
from the optical spectra. As the parameters are interrelated , 
an iterative proces s is used (see, for exam ple, iKilianl 1 19921 ; 
iHunter et ail 120051 ; iThompson et all l2006t ). Studies of gi- 
ants in M 13 yield a cluster metallicity of [Fe/H] = -1.6 
l|Yong. Aoki fc Lambert| [2006). while for uj Cen, which con- 
tains at least three stellar populations, the main cluster g i- 
ant population indicates [Fe/H] = -1.6 (|Origlia et al.l l2003). 
Hence, for both stars the grid corresponding to low metallic- 
ity regions (6.4 dex; -1.1 dex lower than Galactic) was used, 
as this was the lowest metallicity model available and thus 
the most approp riate based on t he metallicities of the re- 
spective clusters. iLee et all (120051 ) found that using grids of 
different iron abundances changed the effective temperature 
estimate by typically ^ 500 K, while the gravity and micro- 
turbulence estimates were unaffected, so the use of this grid 
should not be a major source of error. The adopted atmo- 
spheric parameters and associated errors for both stars are 
listed in Table □ 

Effective temperature estimates, T e g, were determined 
using the silicon ionization balance, with Si n/Si ill employed 
for Barnard 29 and Sim/Siiv for ROA5701. Surface gravi- 
ties, log g, were derived by over-plotting the observed spec- 
trum in the regions of the Balmer lines with theoretical pro- 
files, with associated errors assessed from the uncertainty 
in the fitting (Fig. [l|. The microturbulence, £, was deter- 
mined using the observed O II and Si ill lines, by removing 
the dependance of abundance on line strength, i.e. a plot 
of abundance against line strength having a zero gradient 
l|Kilian|[l992l ). Individual multiplets were employed in this 
procedure, to reduce errors which may a rise from using many 
different multiplets jHunter et alj|2005l ). 

Initial atmospheric param eters for Barnard 29 were 
taken from IConlon et all (| 19941 ), and are T cff = 20000 K, 
log g = 3.0 and £ = 10 km s _1 . The temperature estimate 



remained unchanged at 20000 K with an associated error 
of ±1000 K, the gravity was constrained to log g = 2.95 ± 
0.10, while a determination of the microturbulence led to a 
value of £ = 4 ± 2 km s" 1 . 

For ROA 5701, in itial parameters were taken from 
IThompson et all (|2006l ). namely T cff = 25000 K, log g = 3.25 
and £ = 5 km s _1 . The present analysis led to no change 
in the effective temperature. As no absorption line was de- 
tected for Hell 4541 A, assuming a normal helium abun- 
dance also allowed us to place an independent upper limit 
of T e fi < 26000 K. Therefore, the associated error for T e a is 
±1000 K. The gravity was determined to be log g = 3.3 ± 
0.1, and the microturbulence £ = 1 +j km s _1 . 



3.2.1 Rotational velocity, v sini 

In addition to the microturbulent broadening, some rota- 
tional broadening may occur. The effect of this was investi- 
gated by fitting instrumentally broadened non-LTE theoreti- 
cal profiles for the Si ill multiplet at 4560 A, calculated using 
the appropriate T c ff , log g and £ values, to the observed spec- 
tra. The theoretical profiles were rotationally broadened and 
over lain in order to det ermine the most appropriate vsim 
(see lTrundle et al.|[200ll for further details). For Barnard 29, 
the profile fitting resulted in a value of vsim = km s -1 , 
while in the case of ROA 5701, a vsim = 3 km s~ x was 

obt ained. 

iLehner et ail (2000) discuss the difficulty in differentiat- 
ing between the effects of the two types of broadening when 
such low values are obtained. In order to place an upper 
limit on the vsim values, a microturbulence of km s _1 
was assumed for both stars, as a non- zero microturbulenc e 
would lead to lower vsim estimates (|Trundle et al.ll200ll ). 
This resulted in, for Barnard 29 vsim ^ 6 km s _1 and for 
ROA 5701, vsim sC 3 km s" 1 . 



3.3 Abundances from optical spectra 

Abundances for both stars were determined using the 
adopted atmospheric parameters; atomic data were taken 
from http://star.pst.qub.ac.uk Table [2] contains the equiv- 
alent width measurements for the observed metal lines, 
along with the associated abundance estimates. If more than 
one line contributed to an equivalent width, the lines were 
treated as blends within tlusty. 

There are currently no non-LTE grids available for 
Aim, S ill or Fein, although sulphur is included in our non- 
LTE model atmosphere calculations using a relatively simple 
model. Using tlusty, models were calculated, adopting the 
appropriate atmospheric parameters. These were then em- 
ployed to generate theoretical LTE equivalent widths (see 
IHunter et alj [2005) , which were used to derive the abun- 
dance estimates shown in Table with errors calculated by 
utilising the same method as for other elements. 

Table [3] contains a summary of the abundances 
for each species, along with the associated error esti- 
mates and the number of lines observed. Also included 
are the abundances previously derived f o r the objects 
I Conlon et al. 1 1 1994 iDixon fc Hurwitdll998l ; iMoehler et al.l 



19981 ; Thompson et al. 2006h . and those of n ormal B-type 



stars, both in the Galactic cluster NGC6611 (jHunter et al.l 
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Table 2. Absorption line equivalent width measurements from optical spectra and derived abundances for Barnard 29 and ROA 5701. 



Identification Barnard 29 ROA 5701 

Rest A Species Observed A W Quality^ Abundance Observed A W Quality^ Abundance 
(A) (A) (mA) (A) (mA) 



4267 00 


C n a 




<10 




6578 05 


C n 




<10 




3995 00 

kjUUxJ' KJKJ 


N o a 


3995 00 

kjUUxJ- KJKJ 


65 




403^ 08 


N n 


403^ 0Q 


12 




4041.31 


N II 


4041.32 


25 




404^ ^ 


N II 


4043 ^9 


14 


b 


4227.74 


N II 


4227.71 


g 




4241.78 


N n a 


4241.79 


21 




4432.73 


N o a 


4432.72 


12 




4447 0^ 

'Ttl'rt ( . KJrJ 


N II 


444fi QQ 


30 




4fi01 48 


N II 


4fi01 47 


27 




4607.16 


Nn 


4607.14 


21 


b 


461 3 86 


Nn 


461 3 85 


18 




4621 29 


Nn 


4621 38 


16 






N II 


4(^30 ^fi 


50 




464^ 0Q 


Nn 


4643.10 


27 




^001 in 


N u ab 


^001 19 


30 


b 




N u ab 


5001.47 


37 


b 


^009 70 


N u ab 


^009 fi8 


g 






Nn 


^OO'S 1 ^ 

•J UUu. i-tj 


47 


b 


5007 33 


Nn 


5007 23 

• > KJ KJ 1 . iitj 


30 




3911.96 


O \i ab 








391 2 09 

iW J. — . KJiS 


O n" 6 








3919 28 


O ii 








^4^ 04 


O ii 








QQ^4 37 


O ii 








3Q89 79 


O ii 








40fiQ (S9 


O ii" 6 


4069 61 


13 


b 


40fiQ 8Q 


O ii" 6 


40fiQ 87 


17 


b 


4072.16 


O ii 


4072.13 


23 




407^ 8fi 


O ii 


407^ 83 


30 




407K 84 


O ii 








4089 9Q 


O ii 








41 39 80 


O ii 








41 8^ 4fi 


O ii 








4294.79 


O ii 








4303 84 


O n a 








4317.14 


n a 


4317.12 


10 




4*31 Q 63 

rttj J. tJ • KJ'J 


O ii 


4319.61 


13 




434Q 43 


O ii 


434Q 43 


25 




43M 0(\ 


n a 


43M 9^ 


9 




43fifi QO 


n a 


43fifi 8fi 


12 


b 


43fiQ 97 


O ii 








4oyo.y4 


n tt 

KJ 11 








4414.91 


On 


4414.90 


26 


a 


4416.98 


On 


4416.98 


20 


a 


4452.38 


On 








4590.97 


On 


4590.98 


11 


b 


4596.17 


On 


4596.17 


8 


c 


4609.44 


On a 








4638.86 


On 


4638.87 


14 


a 


4641.81 


On 


4641.83 


27 


a 


4649.14 


Qn ab 


4649.14 


39 


a 


4650.84 


Qn ab 


4650.84 


14 


a 


4661.64 


On 


4661.63 


18 


a 


4673.74 


On 








4676.24 


On 


4676.21 


16 


b 


4699.00 


Qn ab 








4699.22 


Qn ab 








4705.35 


Oil 








4710.01 


On a 









<6.13 ••■ <10 •■• <6.27 

<6.64 ••■ <10 •■• <6.54 



7.22 


3995.00 


45 


a 


6.99 


7.62 










8.30 










7.47 










7.36 










7.62 


4241.82 


15 




7.67 


7.39 


4432.75 


\\ 




7.35 


7.28 


4447.03 


18 


b 


6.95 


7.34 










7.29 










7.34 










7.12 










7.28 




26 


b 


6.93 


7.33 


4643.10 


H 




6.93 


7.36 


^nm ifi 


20 




6.92 




5001 50 

•J KJ KJ _i_ . <J KJ 


25 






7.34 


5005 1 7 


30 




6.88 


7.68 


5007 34 


l\ 


b 


6.96 




391 1 98 


27 




7.64 




3919 30 

J. i7. <JKJ 


19 


b 


7.74 




304^ 04 


15 


b 


7.62 




qq^4 38 


29 




7.65 




3982 72 


20 


b 


7.78 


7.55 


4069 64 


36 




7.76 




4069 90 

t:UU J . UKJ 


47 






7.52 


4072.17 


57 




7.72 


7.54 


4075.87 


68 




7.71 




4078 8^ 


19 


b 


7.79 




408Q 9Q 


33 




8.02 




4132.81 


24 


b 


7.76 




4185.47 


24 




7.67 




4294.77 


\\ 




7.56 




4303 S'S 


15 




7.39 


7.35 


4317.15 


37 




7.78 


7.52 


4319.64 


46 




7.95 


7.59 


4^4Q 4^ 


71 




7.96 


7.29 


4351.27 


38 




7.54 


7 /in 


/i Qfifi on 


Q7 
O J 


a 






4369.30 


10 


c 


7.80 




4395.95 


11 


c 


7.69 


7.53 


4414.91 


70 


a 


7.86 


7.64 


4416.98 


52 


a 


7.89 




4452.37 


16 


b 


7.90 


7.41 


4590.99 


43 


a 


7.66 


7.38 


4596.19 


33 


a 


7.61 




4609.44 


12 


c 


7.36 


7.54 


4638.87 


50 


a 


7.87 


7.52 


4641.83 


79 


a 


7.85 


7.53 


4649.16 


102 


a 


7.86 




4650.86 


48 


a 




7.61 


4661.65 


50 


a 


7.81 




4673.75 


13 


b 


7.90 


7.62 


4676.25 


44 


a 


7.83 




4699.03 


34 


c 


7.63 




4699.24 




b 






4705.37 


37 


a 


7.70 




4710.04 


11 


c 


7.93 
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Table 2. (continued) 



Identification Barnard 29 ROA 5701 

Rest A Species Observed A W Quality^ Abundance* Observed A W Quality^ Abundance* 
(A) (A) (mA) " (A) (mA) 



A son OQ 


(~1 TT 










4oyu.o i 




a 


i . /y 




Cl TT 

\J 11 










4yuo.o4 




c 


7 fin 


A no A £ Q 


U II 










A no A r. Q 


28 


a 


7.91 


4941.07 


On 










4941.10 


12 


c 


7.90 


4943.00 


On" 










4943.01 


21 


a 


7.94 


4481.13 


Mgn ai > 


4481.12 


19 


a 


6.10 


4481.16 


10 


b 


6.15 


4481.33 


Mg ii at 


4481.32 


11 


b 




4481.34 


7 


b 




4128.05 


Sin 


4128.04 


11 


c 


6.45 










4552.62 


Sim 


4552.60 


71 


a 


6.45 


4552.63 


67 


a 


6.18 


4567.82 


Sim 


4567.86 


53 


a 


6.43 


4567.86 


50 


a 


6.18 


4574.76 


Sim 


4574.77 


30 


a 


6.48 


4574.77 


23 


a 


6.16 


5739.74 


Sim 


5739.76 


26 


b 


6.59 


5739.78 


23 


b 


5.96 


4088.85 


Si iv 










4088.89 


33 


a 


6.15 


4116.10 


Si iv 










4116.12 


21 


b 


6.17 


4512.54 


Aim 


4512.58 


7 


b 


5.13 










4529.18 


Aim 


4529.18 


11 


b 


5.15 




<5 




<4.87 


4253.59 


S in 


4253.57 


7 


a 


5.72 


4253.60 


14 


b 


5.62 


4419.60 


Fc in 


4419.60 


3 


c 


6.03 










5156.12 


Fc in 


5156.09 


5 


c 


6.08 










5833.93 


Fc in 


5833.94 


4 


c 


6.21 


5833.99 


<6 


c 


<6.10 



' a indicates that the equivalent width measurement should be accurate to better than ±10 per cent, b to better than ±20 per cent, and 
c to less than ±20 per cent. 

* Logarithmic abundance [M/H] on the scale log[H] = 12.00. 
a Treated as blends in TLUSTY. 

b Equivalent widths of adjacent lines were combined to produce the corresponding abundances shown. 



l200tj ) and in the local field (jKilianl Il99i Il994n . The for- 
mer were determined using the same tlusty code as em- 
ployed in the present paper. Also included are typical in- 
terstellar medium (ISM) abundances for the elem ents con- 
cerned, as listed in lwilms. Allen fc McCravl (l2000ll. The C , 



N and O abundances are taken from 



Cardelli et all (Il996l) 



andlMever. Cardelli fc Sofial (|l997l '); lMever. Jura fc Cardellil 
(1 199sh respectively , with all other elements given by 
ISnow fc Witd lll996lh The Al abundance, how ever, is taken 
as 70 p er cent of the adopted solar abundance of lWilms et al.l 
(|2000l ). as ISM abundances appear to be redu ced by 20 - 30 
per cent with respect to sola r abundances (|Snow fc Wit j 
ll996l ; ISavage fc Sembachlll996h . 

The errors associated with the absolute abundances (see 
Table [3| arise from systematic and ran dom errors (full de- 
tails can be found in iHunter et af1l2005t ). Random errors are 
connected with the data analysis, for example: line fitting 
errors, oscillator strengths and observational uncertainties. 
The random uncertainty is taken as the standard deviation 
in the abundances of a given species, divided by the square 
root of the number of lines observed for that species. Where 
only one line was detected, the standard deviation is that of 
the best observed species, e.g. O II. Systematic errors arise 
from the uncertainties in the atmospheric parameters, and 
are estimated by varying each parameter by its associated 
errors. The total uncertainty is the square root of the sum of 
the squares of the random and systematic errors. We note 
that the scatter of individual abundances for each species 
is reasonably small, with any discrepancies probably arising 



from poorly determined lines and possible errors in atomic 
data. 

A normal He abundance has been assumed throughout 
the analysis of both stars, which was tested by over-plotting 
theoretical spectra on the observations in regions containing 
He I absorption features. There is good agreement to within 
the uncertainties of the atmospheric parameters, suggesting 
this is a valid assumption. When lines are badly blended 
with hydrogen or helium lines, they are not included in the 
non-LTE calculations. 

For both stars, the best observed species are Nil and 
O II, with the majority of equivalent widths believed to be 
accurate to better than 20 per cent. If the lines of poorer 
quality are removed, the abundance estimates remain un- 
changed (within the error bars), indicating that data of 
poorer quality do not significantly affect the abundances ob- 
tained. 

The C abundances for both Barnard 29 and ROA 5701 
in Table [3] are quoted as upper limits. Due to the quality of 
the spectra, it is was only possible to place upper limits on 
the equivale nt widths of the Cn lines at 4267 and 6578 A. 
ISigutl HHl) finds that the 4267 A multiplet gives a limiting 
precision of ~ ± 0.2 dex, over a T e s range of 15000 - 31000 
K (as reproduced bv lNieva fc Przvbilld 120061 ). However, for 
T c ff < 25000 K, the line at 6578 A is somewhat more reliable 
than that at 4267 A. Therefore, the upper limits given in 
Table [3] are based on the 6578 A line alone. 

The Mgn doublet at 4481 A was observed and treated 
as a blend within tlusty, and for both stars the derived 
abundances appear reliable. We note that the Si abundance 
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Tabl e 3. Absolute abundances for Barnard 29 and ROA 5701 fr om t heir optical spec t ra, a long with those previously derived by 
IConlon et al 1 l|l994h : lDixon fc Hurwitj [119981) ; iMoehler et al.l (jl998l l and iThompson et al.l l[ 200y h Also listed are abundances for normal 
B-type stars fro m the local field jKi lian 1992, 1994ft an d the Galactic cluster NGC661 1 (Hunter ct al. 2006J), plus values for the inter- 
stellar medium aCardelli et al,l ll996; Meyer ct al. 1997, 1998; Snow fe Witti ll996; Wilms et al. 2000). Bracketed values are the number 
of lines of each species observed. 



Species Barnard 29 ROA 5701 B-type ISM 

This paper Con 94 Dix98 Moe98* This paper Moe98* Tho06 Kil 92/94 Hun 06 



Cn 


<6.64 


(1) 


<6.70 


6.15 


<6.54 


(1) 


<5.85 


<6.51 


8.20 


7.95 


8.38 


Nil 


7.43 ± 0.12 


(17) 


7.30 




7.06 ± 0.12 


(9) 


6.86 


7.05 


7.69 


7.59 


7.88 


On 


7.50 ± 0.24 


(17) 


7.60 




7.76 ± 0.10 


(41) 


7.97 


7.75 


8.55 


8.55 


8.69 


Mgn 


6.10 ± 0.15 


(1) 


6.00 




6.15 ± 0.21 


(1) 




6.27 


7.38 


7.32 


7.40 


Si II 


6.45 ± 0.20 


(1) 




















Si in 


6.45 ± 0.26 


(3) 


6.30 




6.17 ± 0.20 


(3) 


6.14 


5.94 


7.28 


7.41 


7.27 


Si iv 










6.16 ± 0.40 


(2) 




5.95 








Aim 


5.14 ± 0.10 


(2) 


5.19 




<4.87 


(1) 




<5.21 


6.22 




6.33 


Sin 


5.72 ± 0.26 


(1) 


<6.29 


5.34 


5.62 ± 0.23 


(1) 




5.71 


7.25 




7.09 


Fein 


6.07 ± 0.12 


(3) 


<6.70 


5.28 


5.30 <6.10 


(1) 


4.79 






7.50 


7.43 



* average of curve-of-growth and spectrum synthesis values. 



estimates are sensitive to the adopted atmospheric parame- 
ters. For example, for Barnard 29 a change in T c g of ±1000 
K leads to a variation of ±0.15 dex in the abundance of Sill 
and ±0.21 for Sim. In the case of ROA 5701, the same tem- 
perature change produces an abundance variation of ±0.10 
dex for Sim and ±0.34 dex for Siiv. 

The Sim line at 5739 A was observed in both stars. 
However, as there are no non-LTE grids available for this 
spectral region, theoretical LTE equivalent widths were gen- 
erated, as for S in, Aim and Fein. The resulting LTE 
abundances are 6.59 dex for Barnard 29 and 5.96 dex for 
ROA 5701, which compare w ell t o the S i m re sults in Ta- 
ble [3] [Gie^&Ta^erJ(ll992) and lKilianl (|l994l ) found that 
LTE abundances agree to 0.2 dex with corresponding non- 
LTE abundances for B-type stars, consistent with the results 
here. However, as the 5739 A abundances are in LTE, they 
are not included in the relevant non-LTE Si m results listed 
in Table El 

The S in 4254 A line was used to derive the S abundance 
in both Barnard 29 and ROA 5701. As previously noted, 
these abund ances are based on t heoretical LTE S ill equiv- 
alent widths. iHunter et al.l l|2005l ) find that non-LTE effects 
appear to be negligible for this species. Additionaly, we have 
carried out test calculations and find little difference be- 
tween the S ill spectrum predicted using LTE or non-LTE. 
Hence our use of an LTE approach is probably valid. 

Two Aim lines, at 4512 and 4259 A, were observed for 
Barnard 29. No Al features were detected for ROA 5701, but 
an upper limit on the Al abundance was set using the 4529 
A line. The Al abundances are t reated in LTE, but should 
be reliable as IHunter et all (2005) state that when using the 
dominant ion stage, LTE abundance calculation will give 
good agreement with non-LTE abundances. 

Three Fern lines were observed at 4419, 5156 and 5833 
A in Barnard 29, with one in ROA 5701 at 5833 A. The 
S/N of the local continuum for the lines are ~130, 155, 
150 and 110 per pixel, respectively. LTE calculations were 
used to determine the Fe abundance, and Fig. [2] shows the- 
oretical profiles corresponding to the measured abundances, 



over-plotted onto the observed spectra of Barnard 29 (a,b 
and c) and ROA 5701 (d). As can be seen from the figure, 
some of the features are difficult to distinguish from the 
noise. Indeed, Fig. [U^a) appears to show the 4419 A line 
in Barnard 29, but numerous noise features of similar size 
to the 4419 A line can also be seen, which may be due to 
observational effects or the reduction process. Unless these 
additional features are real (which is unlikely), then the Fe ill 
line may also not be real and hence the abundance derived 
should be used with care. In addition, the 4419 A line ap- 
pears to be narrower than other weak stellar absorption 
features in the spectrum, as may be seen from the profile 
fit in Fig. Ela) , where the theoretical feature is wider than 
that observed. However, 4419 A gives a similar abundance 
to the apparently more reliable 5156 A line, and so all three 
observed Fern transitions are used in determining the iron 
abundance for Barnard 29. Also, these lines have the largest 
gf-values in their respective multiplets, and are among the 
strongest unblended Fe ill lines in the optical spectral region. 
For ROA 5701, the sole Fein line observed at 5833 A does 
not appear reliable (see Fig.[3Jd)) and hence the abundance 
is taken as an upper limit. 

As the observed optical iron lines are of low quality, ul- 
traviolet spectra were also employed to determine Fe abun- 
dances for both Barnard 29 and ROA 5701. This is discussed 
below. 

3.4 Iron abundances from ultraviolet spectra 

To determine an iron abundance from the GHRS/HST 
spectrum, a goodness-of- fit method was used l|Rvans et al.l 
120021 ; iDufton et all I2OO6I , Dufton (in preparation)). Theo- 
retical spectra were generated for the regions observed by 
the GHRS /HST, using non-LTE tlusty model atmospheres 
and appropriate stellar atmospheric parameters. The model 
grid for low metallicity regions was again used, and the un- 
derabundance of elements (heavier than He, up to and in- 
cluding Zn), was varied in each spectrum (in steps of 0.1 dex 
from 0.0 to -2.5 dex) relative to the solar estimates adopted 
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Figure 2. Observed and theoretical line profiles for Fein in Barnard 29 (a,b and c) and ROA5701 (d). The theoretical profiles (smooth 
lines) are the best fits to the observed features, and have been calculated at abundances of (a) On 4414 A = 7.53 dex, On 4416 A = 
7.64 dex, and Fein 4419 A = 6.03 dex; (b) Fein 5156 A = 6.08 dex; (c) Fein 5833 A = 6.21 dex, and (d) Fein 5833 A = 6.10 dex. 



in tlusty |Grevesse fe Sauvalll 19981 ) . Light elements (C, N, 
O, Mg, Si) were treated in non-LTE while all other elements 
were in LTE. Dufton (in preparation) found that ~80 per 
cent of absorption in the ultraviolet spectral region covered 
by our GHRS/HST data are due to Fe, with more than 95 
per cent from Fe-group elements, therefore the calculations 
are taken as being in LTE. Iron dominates the absorption, 
therefore any variation in the relative element abundances 
should not significantly affect the theoretical spectra. 

Due to the numerous absorption features in the ul- 
traviolet spectral region, continuum placement was diffi- 
cult, as may be seen from the GHRS/HST observations of 
Barnard 29 and ROA5701 in Fig. Therefore, the contin- 
uum was not normalized, instead the flux calibrated spec- 
trum was examined. This was achieved by smoothing the 
theoretical spectra with a Gaussian function (to allow for 
instrumental broadening), followed by binning to the exact 
wavelength range and interval as the observed spectra, and 
then scaling to the same mean flux as that observed. 

A goodness-of-fit was calculated by finding the sum of 
the squares of the difference between the observed and the- 
oretical spectra for a range of metallicities. A plot of this 
produced a parabolic shape, with the model with the best 
goodness-of-fit situated at the minimum of the parabola. 
This method removes the need to normalize the observed 
spectrum, thus reducing errors associated with fitting the 
continuum. To test if different spectral intervals would affect 



the abundance estimate, smaller wavelength ranges were se- 
lected and analysed. No significant change was found in the 
goodness-of-fit and associated metallicity value. 

Results of the minimization procedure ar e shown in Ta- 
ble [4] alon g with iron abundances d erived by iMoehler et alj 
(1998|) and bixon fc Hurwitzl l| 19981 ). Errors were calculated 
using a similar method as for the optical spectra. Systematic 
errors were obtained by varying the atmospheric parameters 
by their respective errors, while random errors were esti- 
mated by over-plotting theoretical spectra around the mini- 
mum and observing where the spectrum no longer matched 
the strongest observed features. The total uncertainty was 
then calculated and is shown in Table [4] 

The theoretical spectra that produced a minimum are 
shown, together with the observed spectra, in Fig. [3] For 
Barnard 29, the goodness-of-fit indicates an [Fe/H] = -2.1. 
In general, the theoretical spectrum matched many of the in- 
dividual absorption features in the observational data, sup- 
porting the abundance extimates for this star. ROA 5701 
was less convincing in the fitting of the observed spectrum, 
with a minimum produced for [Fe/H] = -2.4 dex. In par- 
ticular, the spectrum has fewer identifiable absorption fea- 
tures compared to Barnard 29, making the estimation of the 
metallicity unreliable. 

In Fig. there is an emission feature at ~1883 A in 
the spectrum of Barnard 29, which appears real and resem- 
bles part of a P Cygni profile. However, the apparent line 
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does not have a strong associated absorption feature and is 
therefore thought to be an artefact of the reduction process. 
If it is removed from the goodness-of-fit analysis, there is 
no change in the abundance obtained. The feature was not 
observed in ROA 5701. 



4 DISCUSSION 

Barnard 29 and ROA 5701 have been the subject of numer- 
ous studies, several of which have included determinations 
of atmospheric parameters and chemical compositions. We 
discuss these below in the context of the present results. 



4.1 Comparison of atmospheric parameters 

4.1,1 Barnard 29 

Previous studies of Barnard 29 have produced similar atmo- 
sph eric parameters to tho se found in this paper. For exam- 
ple, lAuer fe Norrisl (1 19741 ) used non-LTE models with opti- 
cal spectra to der ive T eff = 2250 ± 1200 K and log g = 3.0 
± 0.15 dex, while Ide Boerl |l985) employed ultraviolet spec- 
tra from the IVE satellite and ultraviolet colours, finding 
T cff = 19000 ± 1000 K and log q = 3.2 dex. T hese spectra 
were again employed by lAdelman et all l| 19941 ) . along with 
optical spectra of H7, to determine LTE atmos pheric param- 
eters of r e ff = 20250 K and log g = 3.15 dex. IConlon et all 
( 1994) used LTE models to analyse IV E low-resolution data 
and WHT/UES optical data, resulting in T cff = 20000 ± 
1000 K, log q = 3 ± 0.1 dex and £ = 10 ± 2 km s _1 . 
iDixon fc Hurwitz 1 19981) too k a sample of the IVE spectra 
used bv lConlon et all (|l99l ). along with ORFEVS-SPAS II 
spectra, a nd with LTE techniqu es obtained the same surface 
gravity as IConlon et all \ 19941 ). but an increased tempera- 
ture estimate of T eB = 21000 ± 1000 K. 

The previously derived effective temperatures and sur- 
face gravities are in good agreement with those found in 
this paper, wi t hin t he errors, with the exception of the 
lAuer fc Norrisl (I1974T ) tempera ture, which may be d ue to us- 
ing different non-LTE models. IConlon et all l| 19941) derived 
a larger microturbulence than in the present paper, presum- 
ably due to u sing LTE models rat her than our non-LTE ap- 
proach. Also IConlon et all (| 19941 ) only employed 11 lines 
in their determination of £, compared to our use of both the 
Sim multiplet and a variety of different multiplets of On. 
Given the high quality of our observational data combined 
with the complexity of our non-LTE models, we believe that 
our atmospheric parameters for Barnard 29 are probably the 
most reliable estimates to date. 



4.1.2 ROA 5701 

ROA 5701 was studied recently by iThompson et all (|2006l ) 
using the same non-LTE models and methods as in the 
present paper, resulting in atmospheric parameters of T e ff 
= 25000 ± 1000 K, log q = 3.2 ± 0.2 dex and £ = 5 ± 5 
km s- MMoehler et all (|l998l ) employed ESO CASPEC op- 
tical spectra, along with the GHRS/HST data used here, 
and found T cB = 22000-24500 K, log q = 3.2-3.4 dex. In 
addition, they derived £ = 2-3 or 20 km s , the former 
using Fein transitions and the latter those of On. Their 



use of LTE models, compared to the non-LTE approach in 
the present paper, may be the reason for the differences 
in the derived microturbulence values. T he present micro- 
turbu lence estimate is smaller than that of IThompson et all 
(2006), probably due to improved equivalent width measure- 
ments and surface gravity estimate, although it is still within 

the uncertainty quoted. 

Other studies of ROA 5701 include INorrisl (|l974l ). who 
analysed narrow-band and UBV p hotometry to obtain T e g = 
26600 K and log q = 3.5 dex, and lCacciari et all l|l984l ) who 
used IVE spectra to derive T c ff = 24000 K. These studies are 
in agreement with the parame ters der i ved h e re, with the ex- 
ceptio n of the T c ff estimate of INorrisl l|l974l ). ICacciari et all 
( 1984.) sugges ted that the different values of (B-V) used by 
INorrisl i| 19741 ) are responsible for their larger temperature 
estimate. The results derived here are hence generally con- 
sistent with previous values, but once again are believed to 
be more reliable due to the improved observational data and 
non-LTE models employed in our analysis. 



4.2 Comparison of chemical compositions 

Table [3] shows the absolute abundances derived for 
Barnard 29 and ROA 5701 , alon g with those previousl y de- 
termined bvlConlon et all (|l994l ); rMo"ehler et all (Il998l) and 
IDixon fc Hurwitj ll 19981) for Barna rd 29, and iMoehler et all 
l| 19981 ) and IThompson et all (|2006h for ROA 5701. However, 
absolute abundances may contain systematic errors, due for 
example to inaccuracies in the atmospheric parameters and 
atomic data, in particular the f-values. In order to mini- 
mize these errors, a detailed line-by-line differential abun- 
dance analysis can be performed, with respect to a young 
B-type star of known chemical composition (Population 1) 
which has similar atmospheric parameters. Therefore, in Ta- 
blc[5]the abundances of Barnard 29 and ROA 5701 are com- 
pared to th ose of young B-typ e stars in the Galactic cluster 
NGC 6611 llHunter et al.ll2006l ) and in the local field (|Kilianl 
1 19921 . 1 19941 ). For the former, results have been determined 
using the same non-LTE model atmosphere calculations as 
in the present paper. The latter uses different model atmo- 
sphere calculations and hence differential abundances us- 
ing this source may not be as reliable. The standard error 
in the differential abundance estimates for Barnard 29 and 
ROA 5701, as compared with young B-type stars, are typi- 
cally ±0.2 dex. These errors do not include any systematic 
errors due to the atmospheric parameters as these are as- 
sumed to cancel due to the use of the same methods for 
both samples. Note that where we did not have access to re- 
sults for individual absorption lines, differential abundance 
values are based on the mean abundance estimates listed, 
with a minimum error of ±0.2 dex. 

It is not possible to perform a differential analysis for 
Fein ultraviolet lines using young B-type stars, as these 
are affected by line crowding and are often saturated. The 
strengths of the Fein absorption lines are only weakly de- 
pendent on the effective temperature and surface gravity, so 
they should not be a significant source of error. In order to 
test the dependence on microturbulence for young B-type 
stars, the Galactic star HR1886 was examined using STIS 
spectra taken from the HST archive. T his star was chosen as 
its T cff = 24000 K and log g = 4.1 dex (|Hamblv et allll996l ) 
are similar to those of ROA 5701. For a microturbulence of 
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Figure 3. The upper panels show the observational data for Barnard 29 and ROA 5701 obtained with GHRS/HST, 
panels contain the theoretical spectra for [Fe/H] = -2.1 and -2.4, respectively, relative to the Galactic value of [Fe/H] = 
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= 7.5. 



Table 4. Values of [Fc/H] derived for Barnard 29 and ROA 5701 relative to the Galactic abundance, along with previous estimates. 



Star 



This paper Moehler et al. (19981 Dixon fc Hurwitz (1998) 



Barnard 29 -2.1 ± 0.1 
ROA 5701 -2.4 ± 0.3 



-2.2 
-2.7 



-2.2 



km 



synthetic spectra were generated and scaled to 



the same mean flux as the observed spectrum, using the 
methods described in Section \3. 41 The microturbulence was 
then increased to 10 km s _1 , and the process repeated. The 
change in microturbulence resulted in a change of metal- 
licity of ~1.4 dex, implying that Fein ultraviolet lines in 
Galactic B-type stars are indeed very sensitive to changes in 
microturbulence, and are not appropriate for a differential 
abundance analysis. 



iKendall et ail (|l994al ) found a similar dependence 
on microturbulence for the standard star HD 149438 
(rScorpii), with a change of £ from to 10 km s yield- 
ing a metallicity change of ~0.8 dex. For Barnard 29 and 
ROA 5701, a £ change from to 10 km s _1 resulted in a 
metallicity change of ~0.4 and 0.2 dex, respectively. 



4.2.1 Barnard 29 



Differential abundances for Barnard 29, relative to young 
B-type stars, are given in Table \5j along w i th th e re- 
sults of IConlon et al l i| 19941 ); iDixon fc Hurwit3 il 19981) and 
iMoehler et al.l 1 19981 ). As stated previously, Conlon et al.l 
( 19941 ) used IV E and WHT/UES data with LTE models in 
their analysis. Their upper limit for the C abundance agrees 
with that found in this paper, while N, O, Mg, Si and Al are 
~0.1 dex lower. Reasons for the differences in abundance 
could include different atmospheric parameters, improved 
equivalent width measurements (as a consequence of our 
higher S/N observations), and our use of non-LTE model 
atmosphere codes. To investigate if our non-LTE treatment 
may be responsible for the differences, we have calculated 
LTE abundances for the elements, with the exception of 
Al which is already in LTE. This results in abundances of 
<5.93 for C, 7.40 for N, 7.75 for O, 5.93 for Mg, 6.78 dex 
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Table 5. Differential abundance s for Barnard 29 a nd ROA5701, relative to those in young B-type stars from iHunt er et al, (2006) 
wher e available, supplemented bvlKilianl ( 199 3.11994^. Also liste d are differential abundances relative to the ISM values jCardelli et al.l 
Il996l: iMever et alJ 19971 . Il998l : ISnow fe Wittl Il996l : IWilms et al.lEoOoh . and the metallicities of M 13* l|Yong et al.ll2006h and uj Cen* 
llOrielia et al.ll20031h 



Species Young B-type stars ISM 

Barnard 29 ROA5701 Barnard 29 ROA5701 

This paper Con 94 Dix98 Moe98 This paper Moe98 Tho06 This paper 



c 


<-1.3 


<-1.3 


-1.8 


<-1.4 


<-2.1 


<-1.4 


<-1.7 


<-1.8 


N 


-0.3 


-0.3 




-0.6 


-0.7 


-0.6 


-0.5 


-0.8 





-1.1 


-1.0 




-0.8 


-0.6 


-0.8 


-1.2 


-0.9 


Mg 


-1.2 


-1.3 




-1.2 




-1.1 


-1.3 


-1.3 


Si 


-1.0 


-1.1 




-1.2 


-1.3 


-1.5 


-0.8 


-1.1 


Alt 


-1.1 


-1.0 




<-1.4 




<-1.0 


-1.2 


<-1.5 


s t 


-1.5 


<-1.0 


-1.9 


-1.6 




-1.5 


-1.4 


-1.5 


Fe (Optical) 


-1.4 


<-0.8 




<-1.4 






-1.4 


<-1.4 


Fe (UV) 


-2.1 




-2.2 


-2.2 -2.4 


-2.7 




-2.1 


-2.1 


Cluster [Fe/H] 




-1.6° 






-1.6* 









Note: Con 94 = IConlon et al ] ||1994 ). Dix98 = iDixon fc Hurwit j l|l99Sft . Moe98 = iMoehler et al.l jl998ft . Tho 06 = ThoniDson et al.l 
||2006| ). 



for Si. The N , O an d Si abundances are still larger than the 
IConlon et al.l (Jl994) abundances, while C and Mg are lower. 

The lower microturbulence used in the present analysis 
may be a further source of the differences. Adopting a mi- 
croturbulence of 10 km s _1 from lConlon et al.1 (ll994T ), LTE 
abundances of 7.30 for N, 7.62 for O and 6.47 dex for Si 
are obtained, while the other elements abundances do not 
signi ficantly change. The values are now in good agreement 
with lConlon et ail |l994), with the exception of C, suggest- 
ing that the microturbulence is a factor. The discrepancy 
for C may be due to the use of different lines. However, the 
results in the present paper should be more reliable due to 
our use of non-LTE techniques and improved observational 
data. 

IDixon fe Hurwit3 (fl998j) used a far-ultraviolet spec- 
trum from ORFEUS-SPAS II, fitting LTE synthetic models 
to the data, similar to the method used here. They deter- 
mined abundances for the el ements He, N, O, A l and S, 
finding results consistent with lConlon et all l| 19941) but with 
larger error bars, suggesting the use of lo wer resolution data 
as the reason for this. The C abundance of lDixon fe Hurwit j 
(1998) was found using two Cm lines, and is -0.5 dex lower 
than the upper limit found here. Their S abundance was de- 
rived from lines of S ill and Siv, with some scatter due to 
blending with nearby features. These authors also observed 
S II lines, but these were not included in the analysis as they 
appeared to be affected by non-LTE. The S abundance of 
Dixon & Hurwitz is ~0.4 dex lower than our value, which 
may be due to the use of different spectral wavelength re- 
gions. 

An iron abundance of [Fe/ H] = -2.2 was determine d 
by bothlDixon fe Hurwit3 l| 19981 ) and lMoehler et all j|l998l ). 
IDixon fe Hurwit j ( 1998h Found it difficult to determine an 
Fe abundance, as the spectrum was contaminated by inter- 
stellar lines and many of the expected Fe transitions were in 
regions of poorly-fitted continuum. They fitted a band of Fe 
lines rather than individual features, similar to the method 
used in this paper for the UV observations, noting that their 



abundance was more than 0.5 dex below the model atmo- 
sphere used to derive the synthetic spectra (i.e. [M/H] = 
-1.5), with tests showing that the mod e l choi ce had little 
affect on the abundance. IMoehler et alJ (1 19981 ) determined 
an Fe abundance using the same GHRS/HST spectrum em- 
ployed here. They a dopted the atmospheric parameters of 
IConlon et al. (1994), and fitted a global continuum to the 
normalized data, noting that an uncertainty of 5 per cent in 
the continuum definition of their spectrum would lead to an 
error of 0.2 dex in the derived abundance. 

In the present paper, an Fein abundance of [Fe/H] = 
-2.1 is obtained from the UV spectrum and [Fe/H] = -1.4 
from the optical. The UV value i s well below the m etallic- 
ity of the cluster ([Fe/H] = -1.6, lYong et al.ll2006l ). but is 
in slightly better agreement compared to the previous UV 
studies. However, our optical result is consistent with the 
cluster metallicity, to within the errors. This discrepancy is 
discussed further in Section [ 



4.2.2 ROA5701 

For ROA5701, there is g ood agreement between the present 
abundances and those of lThompson et all (|2006l ). with C, N 
and O giving the same differential values. The Si abundance 
estimates differ by 0.2 dex, and this discrepancy appears to 
be due to differences in the atmospheric parameters. This 
was tested by calculating non-LTE abundances using the 
equivalent wid ths from this paper an d the atmospheric pa- 
rameters from lThompson et al.1 l|2006l ). The resulting values 
are 5.93 dex for Si ill and 5.88 dex for Siiv, suggesting the 
changed atmospheric parameters are a cause for the differ- 



The upper limit to the C abundance of IMoehler et al] 

(1998) is lower than that derived in this paper, despite us- 
ing the same upper limit value for the equivalent width. This 
difference co uld be due to large non-LTE effects in C II, as 
proposed by ISigutJ l|l996l ) and hence, the upper limit de- 
rived here may be more reliable. Although the N and O 
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Table 6. Absolute abundances for Barnard 29 and ROA5701 
from their optical s pectra. Also tabulated ar e the abundances 
of th e M 13 giants llCohen fc Melendej |2005| ° and lYone et al] 
l2006t) and those of uj Ce n, using the main cl uster giant population 
jSmith et al]|2000l '' and lOrielia et alj|2003h . 



their r espective clusters, namely M 13 llCohen fc Melendej 



2005: 



. .Yong et all 120061 ) and uj Cen (ISmith et al.l ,_ 
Origlia et al .1120031 ) . In general the abundances compare well, 



2000; 



Species 


M13 




w Cen 






Barnard 29 


Giants 


ROA5701 


Giants 


C 


<6.64 




<6.54 


6.86 


N 


7.45 




7.06 




O 


7.50 


7.53 


7.76 


7.45 


Mg 


6.10 


6.14 


6.15 


6.26 


Si 


6.45 


6.38 


6.17 


6.26 


Al 


5.14 


5.75 a 


<4.87 


4.87 6 


S 


5.72 




5.62 




Fc 


6.07 


5.90 


<6.02 


5.92 



abundances of iMoehler et al.l (|l99of ) differ by ~0.2 dex, and 
the Si abundance by ~0.1 dex, from those presented here, 
again this may be due to non-LTE effects. To test this, LTE 
abundances have been derived, using the tlusty grid, re- 
sulting in <6.23 for C, 7.99 for O, 7.10 dex for N and 6.50 
dex for Si. The O abu ndance is in very good agreement with 
IMoehler et all <|l998l ). but the C, N and Si values are still 
higher. The changes could be due to the use of different 
equivalent width measurements and atmospheric parame- 
ters. However, the present results should be more reliable 
due to our non-LTE analysis. 

No abundances for M g, Al or S wer e deri ved by 
IMoehler et all i|l998h . while [Thompson et al.1 < |2006f ) deter- 
mined values for Mg and S, which differ from those found in 
this paper by ~0.1 dex. This is probably due to the improved 
S/N observations, which allow better equivalent width mea- 
surements; changing the atmospheric parameters had little 
affect on these abundances. The listed Al abundance is still 
an upper limit, but due to the improved S/N of the spec- 
trum it was possible to place a revised lower limit on the 
equivalent width. In general, the abundances quoted in the 
present paper should be more accurate due to the improved 
qua lity of the observational data. 

iThompson et all (|2006l ) observed two unidentified fea- 
tures at 4512.41 and 4352.22 A in their data, which were 
obtained using UCLES at the AAT. However, no features 
are observed at these wavelengths in the FEROS spectrum, 
which is of higher resolution and S/N, indicating that the 
line s in the UCLES data may not be real. 

IMoehler et a l. ( 1998) obtained a differential Fe in abun- 
dance of [Fe/H] = -2.7 from the GHRS/HST ultraviolet 
spectrum, compared to [Fe/H] — -2.4 derived here. Both val- 
ues are lower than the met allicity of uj Cen, namely [Fe/H] 
= -1.6 (|Origl ia et al . 2003), although our result is in bet- 
ter agreement with the metallicity of the cluster. An optical 
value [Fe/H] < -1.4 is obtained in the present paper. The 
[Fe/H] abundance of ROA5701 is discussed further in Sec- 
tion 031 



4.3 Comparison with Cluster Giants 

The absolute abundances of Barnard 29 and ROA5701 
are given in Table [6] along with those of giants from 



to within the errors. A notable exception is the Al abun- 
dance of Barnard 29 , which is lo wer than that of the g iants. 
However, lYong et all (|2006l ) and ISneden et all (|2004i ) have 
both shown that the aluminium abundances in the giants 
range from [Al/Fe] ~ - 1.2 dex. The giants with [Al/Fe] 
~ 0.3 have similar O abundances to those found here, sug- 
gesting that the Barnard 29 abundance compares well with 
a selection of the giant population. 



4.4 Evolutionary discussion 

4-4-1 Dredge-up 

IThompson et all l|2006] ) previously studied ROA5701, sug- 
gesting that the observed abundance patterns implied that 
the star had left t he AGB prior to th e onset of the third 
dredge- up, as have lConlon et all |l994) to similarly explain 
the abundance patterns in Barnard 29. Dredge-up occurs in 
stars during their AGB evolution, as products of hydro- 
gen and helium burning are moved from the core to the 
surfa ce, altering the surface composition (|lben fc Renzinil 
1983). Briefly, during the first dredge-up, the surface O sees 
little change while C decreases as N increases. The second 
dredge-up sees C and O mostly converted into N. At the 
third dredge-up, triple-a products are brought to the sur- 
face and C is enhanced relative to O and N. 

Both Barnard 29 and ROA 5701 exhibit abundance pat- 
terns compatible with no third dredge-up. The N abun- 
dances are strongly enhanced relative to both C and O, es- 
pecially for Barnard 29, while O is larger than C. This is 
consistent with the previous studies o f both stars, w i th the 
exception of the lower N abundance of lMoehier et all (|l998l) 
for ROA 5701. Thus the abundances of CNO appear com- 
patible with the suggestion that the third dredge-up may 
not have occurred. 

The low C abu ndance could be ex plained in terms of 
hot bottom burning. iDeroo et al.l l|2005l ) found V435 Oph to 
have a low metallicity and no enhancement of C, suggest- 
ing hot bottom burning as an explanation o f the observed 
abun dances. In hot bottom burning (HBB; iMoonev et al.l 
l200ll ). c is converted to N at the base of the convec- 
tive envelope, wh en temp e rature s become large enough for 
CNO processing. iHerwid (|2005l ) indicates that HBB may 
have an important role at ex t remel y low metallicities. How- 
ever, [D^Antona^&^^azzitelh| (|l996l ) suggest that Population 
I stars of initial mass M ^ 5 Mq, and Population II stars 
of M 3.5 Mq experience HBB. As the initial mass of the 
stars is unknown, no conclusions can be made regarding this 
process. 

4-4-2 Gas-dust separation 

Although the lack of a third dredge-up appears a likely ex- 
planation for the abundance patter ns of the light elements 
in bot h stars, IMoehler et al. | l| 19981 ) and [5 ixon fc Hurwitd 
(l998j) explain the abundances of the heavier elements in 
terms of the gas-dust separation process. The basic premise 
involves grain formation that occurs in the circumstel- 
lar shells of stars during the AGB evolution. Dust parti- 
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cles are removed by radiation pre ssure and the 'cleaned ' 
gas is re-accreted o nto the surface (Mathis fe Lamersl[l99i 
IWaters et allll992ft . This results in abundances similar to 
the ISM, i.e. CNO and S near solar, while other elements 
such as Mg, Al, Si and Ca follow the low Fe abundance 
(Van Winckel et al.l[l995h and see for example the RV Taur i 
stars UYCMa, HP Lyr and BZ Set (iGiridhar et al.ll2005h . 
IParthasarathv. Garcia Lario fc Pottaschl (1 19921 ) found that, 
in the atmospheres of post-AGB stars, metals with a high 
condensation temperature, e.g. Fe, are more likely to con- 
dense into dust grains, compared to elements such as C, N, 
O and S with low condensation temperatures which stay in 
the gas phase with no sig nificant depletions observed. 

iMoehler et al. I lll99Sft suggested that this process oc- 
curred as they fo und a low Fe value, with N and O virtually 
undepleted, while IPixon fc Hurwitj (l998l ) also cite the low 
Fe abundance as evidence for the process. For this scenario, 
similar depletions of Mg, Si and Fe would therefore be ex- 
pected. Examining the differential ISM abundances (Table 
[5J it is possible to determine if the stars follow the abun- 
dance trends seen in gas-dust separation. For Barnard 29, 
we find differential abundances of [Mg/H] = -1.3, [Al/H] = 
-1.2, [Si/H] = -0.8, [S/H] = -1.4 and [Fe/H] (optical) = -1.4 
(relative to the ISM). In addition, for ROA5701 the differ- 
ential abundances are [Mg/H] = -1.3, [Si/H] = -1.1 and 
[S/H] = -1.5. The Mg, Al and Fe (optical) abundances of 
Barnard 29 agree with the value for S and the metallicity an d 
giant abundances of M 13 ([Fe/H] = -1.6. lYong et alj|2006tl . 
although the Si abundance estimate is somewhat larger, pos- 
sibly due to its sensitivity to the atmospheric parameters. 
For ROA5701, the S abundance agrees (to within the er- 
rors) with the Mg and Si values and is consis tent with the 
meta llicity and giants of o;cen ([Fe/H] = -1.6. lOriglia et all 
2003). Similar results would be obtained if the abundances 
were compared with young B-type stars instead of the ISM. 

For gas-dust separation, CNO and S would be expected 
to be relatively undepleted and not in agreement with Mg, 
Si, Al or Fe. However, we find that the S abundance is 
generally in agreement with the Mg, Si, Al and Fe(optical) 
abundances, suggesting that gas-dust separation has not oc- 
curred. Also this process is generally confined to binary sys- 
tems. The possible binarity of our two stars was investi- 
gated using radial velocity measurements. However a liter- 
ature search produced sparse results, indicating no conclu- 
sions can be made. Additionally, the position of the stars on 
log g/ T eS or logL/Teff diagrams l|Conlon et al.lll994l ; iBlockerl 
1 19951 ) suggest core masses for both stars o f ~ 0.55 Mq, using 
log(L/L Q ) = 3.2 ilLandsman et al.ll 19921) fo r ROA5701 and 
log(L/L n ) = 3.23 (Adelman et al.lll994l ) for Barnard 29. 
IConlon et all (|l993D found that stars of this average mass do 
not show an infrared excess, which is supported by the lack 
of IRAS detection for both stars, and no emission features 
observed in Ha for ROA 5701. These results suggest there is 
no evidence of a circumstellar dust shell present for either 
star. Hence the gas-dust separation process appears unlikely 
to be responsible for the observed abundance patterns. 



4.4.3 UV lines 

In this paper, we find a difference of 0.7 dex and <1.0 
dex for Barnard 29 and ROA 5701 respectively, between 
the optical and ultraviolet derived iron abundances. Anal- 



yses of the Fein lines in other hot post-AGB stars, such 
as BD+33° 2642 jNapiwotzki. Heber fc Koppenlll994l ) and 



HD 177566 (Kendall et al.lll994al |b h, show similar Fe abun- 
dance trends to those found hereT lNapiwotzki et al.l (1994T ) 
analysed optical and IUE spectra, and their results led them 
to suggest the occurrence of gas-dust separation, as He was 
near solar, C, N, O, Mg and Si d e pleted by ~1.0 dex, and Fe 
depleted by ~2.0. lKendall et all (|l994bl ) set upper limits for 
the Fe abundance from Fe ill lines in both optical and ultra- 
violet spectra, finding a di fference of ~0.5 dex b etween the 
two sets of data. Similarly, IMoehler et al.l (l998ft derived an 
Fe abundance of [Fe/H] = 6.95 dex from the Fe ill lines in the 
IUE spectra of the main-sequence B-type star 7 Peg, com- 
pared to [F e/H] = 7.56 dex from optical Fe hi transitions. 
In addition. iGrigsbv. Mulliss fc Baen (|l996l ) in their analysis 
of iHer found that the ultraviolet Fe abundance was ^0.5 
dex below the near so lar v alue derived from optical lines. 
iGies fc Lambert! <| 19921 ) and lProifitt fc Quiglevl pOOlh stud- 
ied the B-type supergiant HD 51309. While the former state 
two different iron abundance estimates depending on the 
temperature structure used, there is a difference of ^0.75 
dex found between their optical results ( from Fe 11 and Fe hi 
lines) and the IUE iron abundance of IProffitt fc Quiglevl 
(2001ft . 

Dufton (in preparation) have observed similar abun- 
dance discrepancies for young B-type stars in the Small and 
Large Magellanic Clouds (SMC, LMC), as well as in the 
Magellanic Bridge (MB). They determined Fe abundances 
from GHRS/HST and STIS/HST spectra, using the same 
methods and Fe ill transitions as described here, and found 
that these were consistently lower (more than 0.5 dex) than 
the estimated metallicities of the SMC, LMC and MB re- 
gions. 

It is possible to derive the abundance of other elements 
using ultraviolet spectra, but that is not possible here as the 
GHRS data employed do not cover a large enough spectral 
range. However, IUE spectra have been used by a number of 
authors to derive th e abundances of several light elements. 
IConlon et al.l (l99ll ) obtained a C abundance using both ul- 
traviolet an d optical spectra, find ing the abundances to be 
compatible. iKendall et all (|l994bl ) found similar differential 
abundances for C, N and Si between their optical and IUE 
spectra, with differences b eing less than 0.15 dex. A lso, the 
Mg abundance obtained bv lNapiwotzki et al.l (1 1994ft showed 
a difference of 0.05 dex between their optical and ultravio- 
let values. Howev e r, the ultraviolet abundances obtained by 
iNapiwotzki et al.l (1994T ) for C and S are lower by 0.5 and 0.4 
dex, respectively, than the optical abundances determined 
here (see Section ^. 2. 1[) . Despite this result, the use of ultra- 
violet spectra to determine the abundances of light elements 
generally appears to be reliable and do not suffer the same 
problem observed with iron. 



4-4-4 Oscillator strengths 

The above studies, in conjunction with the present work, 
strongly indicate that the ultraviolet Fein lines yield 
systematically low Fe abundances, possibly as a re- 
sult of errors in the adopted atomic data, in particu- 
lar f- values. Currently, ato mic data from Kurucz models 
(http: / /nova.astro. umd.e du , lHubenvlll988l ; iHubenv fc Lanzl 
1995; Hubenv et al.|[l998l ) are used in the tlusty grids. We 
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have therefore performed a literature search of f-values for 
Fein, to investigate if more recent calculations or measure- 
ments would produce different Fe abundances. Table [7] lists 
the six strongest transitions observed in the GHRS/HST 
spectra, along with the oscil lator strengths from Ku rucz, 
and the more recent res ults of lNahar fe Pradhanl l|l996l ) and 
llbner fe HibbertJ (|2005l ). There are little differences between 
the sets of atomic data, so these will not solve the Fe abun- 
dance problem. However we note that further calculations 
and measurements for these Fe ill transitions are planned in 
the future. 

Despite numerous studies of the ultraviolet Fein tran- 
sitions, the line list for the spectral region employed in this 
paper may not be complete. However, it is assumed that, if 
this is the case, any lines excluded during the fitting pro- 
cedure would be compensated by making other transitions 
too strong. This would potentially lead to an over estimate 
in the derived abundance, which would not solve the abun- 
dance discrepancy observed. 

4-4-5 Further work 

To investigate the Fe abundance discrepancy further will 
require very high S/N data for Barnard 29 and ROA5701, 
sufficient to ascertain the reliability of the observed opti- 
cal Fein lines and observe further Fein lines in both stars. 
Unlike the ultraviolet transitions, the optical lines are suf- 
ficiently weak, even in main-sequence Population I B-type 
stars, for a detailed line-by-line differential abundance anal- 
ysis to be performed. This will hence allow the removal of 
systematic effects such as the dependence on f-value, and 
hence the determination of accurate Fe abundances, for com- 
parison with the ultraviolet measurements. We plan such 
observations in the near future. 



5 CONCLUSIONS 

The abundance patterns observed for both Barnard 29 and 
ROA5701 indicate that the stars have not undergone the 
third dredge-up. In addition, the gas-dust separation pro- 
cess suggested by previous authors appears unlikely. The 
optical spectra imply abundances for the Fe and other ele- 
ments, which are compatible with those found in the late- 
type giants, while the Fein lines in the ultraviolet produce 
lower abundances than expected, as found in several other 
studies. This strongly suggests that the ultraviolet Fein do 
not provide reliable abundance indicators. Future studies re- 
quire high S/N optical data, in order to obtain accurate Fe 
abundance estimates. A re-evaluation of the Fein f-values, 
either theoretically or experimentally, would also be highly 
desirable. 
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